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ABSTRACT 

Previous studies of the active galactic nuclei (AGN) contribution to the cosmic 
X-ray background (CXB) consider only observable parameters such as luminosity 
and absorbing column. Here, for the first time, we extend the study of the CXB 
to physical parameters including the Eddington ratio of the sources and the black 
hole mass. In order to calculate the contribution to the CXB of AGN accreting 
at various Eddington ratios, an evolving Eddington ratio space density model is 
calculated. In particular, Compton thick (CT) AGN are modeled as accreting at 
specific, physically motivated Eddington ratios instead of as a simple extension 
of the Compton thin type 2 AGN population. Comparing against the observed 
CT AGN space densities and logA-logS 1 relation indicates that CT AGN are 
likely a composite population of AGN made up of sources accreting either at > 
90% or < 1% of their Eddington rate. 

Subject headings: galaxies: active — galaxies: quasars: general — X-rays: diffuse 
background 



1. Introduction 

It is believed that all forms of active galactic nuclei (AGN) ar e variations of the same 



phenomenon, namely accreting supermassive black holes ( lReeslll984j ). The most basic obser- 
vational description of AGN is their bolometric luminosity, Lboi, and the obscuring column 
density, Ah, through which they are observed. AGN are known to radiate at luminosities 
Lt, i « 10 40 -10 48 erg s" 1 , but this is more physically described by the Eddington ratio, de- 
fined as Lboi/LEdd, where Lem = ^nGM m m p c/ o~t is the Eddington luminosity of a black 
hole with mass M m . Observed Eddington ratios range roughly ten orders of magnitude with 
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a maximum of ~1 fjCao fc Xul 120071 ). There is a similarly wide range of observed column 
densities (20 < logN H < 25; e.g., Tueller et al. 2008). 

The most widely accepted theory to explain the observed range of obscuration is known 
as the unified model and claims that the obscuration of an AGN is dependent solely on the 
orientation of the AGN with respect to the line of sight to the observer (e.g., Antonucci 
1993). However, there is evidence that suggests that Compton thick (CT) obscuration, 
when logiVn > 24, is linked to an evolutionary phase during which the black hole and the 
host stellar spheroid accrete rapidly (e.g., Fabian 1999; Page et al. 2004; Rigopoulou et al. 
2009). As the black hole grows, the outflow from the black hole strengthens and ejects the 
obscuring material, revealing an unobscured quasar. This outflow also halts the growth of 
the host spheroid causing the observed proportionality between black hole mass and host 
spheroid mass (e.g., Crenshaw et al. 2003; Page et al. 2004; Alexander et al. 2010). Galaxy 
merger simulations show that mergers between gas rich galaxies will cause gas to flow into 
the nuclear region, igniting a nuclear star burst and high l y obs cured, high Eddington ratio 
quasar activity (e.g., Hopkins et al. 2006). iFabian et al. I ( 120091 ) also show that if an AGN is 
accreting at close to Eddington, radiation pressure on dust will blow out any column density 
with logiVn < 24. Therefore, if an AGN is obscured and has a high Edding ton ratio, it 



must be CT. Observational evidence also supports this evolutionary scenario. iPage et al. 



( 120041 ) found that star formation in the hosts of unobscured quasars has already peaked, 
whereas the hosts of obscured quasars are still undergoing massive amoun t s of sta r formation 



indicating that their galactic spheroids are still forming. iTreister et al. I (|2009aJ) found that 
the space density of luminous CT AGN candidates evolves strongly at z = 1.5-2.5, so that 
the CT quasar population seems to peak at a slightly higher redshift than the population 
of unobscured quasars. This evolutionary theory is difficult to observationally test since the 
high levels of obscuration found in CT AGN make them nearly invisible below 10 keV in 
their rest-frame. 

If being CT is an early evolutionary phase of powerful quasars, CT AGN should be 
very rare at z ~ 0; however, Swift/BAT and INTEGRAL have detected several CT AGN at 
z < 0.025 (Tueller et al. 20 08; Malizia et al. 2010 ). According to black hole masses and X- 
ray luminosities reported by lGultekin et al. I (120091 ) . the two brightest CT AGN, the Circinus 
Galaxy and NGC 4945, have A = —2.0 and A = —5.6, respectively, where A = log L bo i/L Edd . 
Therefore the scenario where CT AGN are high Eddington ratio quasars does not completely 
describe the entire CT AGN population. 

As the majority of the cosmic X-ray background (CXB) < 10 keV has been resolved 
into AGN by deep observations by ROSAT, Chandra, and XMM- Newton (see Brandt & 
Hasinger 2005), the CXB provides a complete census of AGN. Due to the observational 
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constraints of highly obscured sources, CT AGN are difficult to study even in the local 
universe; thus, historically CT AGN have been invoked to augment the observed AGN 
population in population synthesis models to make these models properly fit the observed, 
but not resolved, peak of the CXB around 30 keV. As there is a dearth of observational 
information about CT AGN, previous population synthesis models have treated CT AGN as 
a simple extension of type 2 AGN (e.g., Ueda et al. 2003; Treister & Urry 2005; Ballantyne 
et al. 2006; Gilli et al. 2007); however, if CT AGN are part of an evolutionary sequence, 
then they would not evolve in the same manner as the less extremely obscured type 2 AGN. 
Indeed, fitting the CXB in that manner seems to overpredict the local observed space density 
of CT AGN (see Treister et al. 2009a). 

In this letter, we compute the contributions of AGN with various Eddington ratios to 
the CXB and model CT AGN using a physically motivated Eddington ratio distribution. We 
then compare our model predictions with observed cumulative and decadal space densities 
and the Swift/3 AT CT logAMogS" relation. We find that CT AGN are a composite popu- 
lation of AGN accreting at both very high Eddington ratios and very low Eddington ratios. 
When necessary, a ACDM cosmology is assumed with h = 0.7, Qm = 0.3, and = 0.7. 



2. Calculations 



2.1. Eddington Ratio Distribution 



Following the method of iMerlonil (120091 ) and iMerloni fc Heinzl ( 120081 ). with minor modi- 
fications described below, we calculate the Eddington ratio space density, (Lx, M„ z), in 
Mpc~ 3 , where Lx is the 2-10 keV luminosity. The bolometric correction from lMarconi et al. 
(b004h . 

\og(L bo i/L x ) = 1.54 + 0.24C + 0.012£ 2 - 0.0015£ 3 , (1) 
where £ = \oz{L hnl /L^ 



nosities. 



Marconi et al. 



— 12, is used to transform X-ray luminosities to bolometric lumi- 
( 120041 ) derive this bolometric correction using an average intrinsic 
spectral energy distribution template for radio quiet AGN. The scatter in the bolometric cor- 
rection is ~0.1 dex, which is negligible in our analysis. The black hole mass function (BHMF) 
is derived f rom a Gaussian curve fit to the combined type 1 and type 2 AGN BHMF observed 
bv lNetzerl J2009h at z ~ 0.15 . As the iNetzerl J2009h BHMF is fractional, the IMerlonil J2009h 
and IMerloni &: Heinzl (120081 ) method was modified to use only the space density from the 



Ueda et al. I (120031 ) hard X-ray luminosity function (HXLF), instead of dep ending equally on 
the space density from the BHMF and the HXLF. IMerloni fe Heinzl (120081 ) argue that black 
hole mergers will only have a perturbative effect on black hole growth as the uncertainties 
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in the AGN bolometric luminosity function are larger than the impact of black hole mergers 
on the local BHMF. Therefore, the evolution of the BHMF must satisfy the conservation 
equation (Small & Blandford 1992): 



dn M (M„ t) + d[n M (M.,t)(M(M.,t))] = Q 



Of 



DM 



(2) 



where tim is the BHMF and (M(M,,t)) is the average accretion rate at time t. Thus the 
BHMF is evolved by integrating this conservation equation forward and backward in time, 
using the observed BHMF at z 0.15 as the boundary condition. The critical Eddington 
ratio where the transition from radiatively efficient thin dis k accretion to ra diatively inef- 
ficient geometrically thick accretion is taken as A cr = —2.0 (ICao fc Xul 120071 ). For A > A cr 



the radiative efficiency is assumed to be e ra d = 0.1 and for A < A cr , e ra d = 0.1(A/A cr ) 1 ' 2 



(Merloni fc Heind 120081 ). 



2.2. Cosmic X-ray Background Model 



The CXB model closely follows that described by iDraper fc Ballantynd (120091 ) with 
differences in the spectra of CT sources and low Eddington ratio sources, and §\(L X , M, z) 
is input i n place of a luminosi ty function. A torus reflection component, calculated using 
"reflion" jRoss fc FabiarJ 12005^ . is added to the CT AGN spectrum (cf., Gilli et al. 2007). 
For AGN with A < —2.0, the reflection component is included only if the source is CT. In 
order to properly fit the CXB at higher energies, the power-law component of the A < —2.0 
spectrum is assumed to have a cut-off energy of 200 keV, otherwise a cut-off energy of 250 
keV is used3. Most importantly, CT AGN are added at specific Eddington ratios which are 
chosen based on a physically motivated scenario instead of adjusting the normalization of 
the luminosity function to account for CT AGN as done with previous CXB models. 



3. Results 

In previous population synthesis models the number of CT AGN was assumed to be 
proportional to the number of type 2 AGN. In this original model we find 44% of all type 2 
AGN are CT, where the normalization is fixed by fitting the peak of the CXB (EFe ~ 44.2 
keV cm -2 s _1 str -1 at 30 keV). In figure [1] the cumulative space density is shown for CT 



1 If 250 keV is used as the cut-off energy for all sources the qualitative findings of this study do not change, 
but the fraction of CT AGN is reduced by ~1%. 
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AGN with logLx > 43 (red) and \ogLx > 44 (blue). The dashed line shows the predictions 
of the original model . The red circ l e show s the local Swift/BAT and INTEGRAL data 



originally reported by iTreister et al. I ( j2009bl ) and here has been c orrec ted to account for the 



flux-luminosity relation for CT AGN described by iRigby et al. I (120091 ) . The original model 
significantly overpredicts the observed space density of CT AGN with logLx > 43 at all 
redshifts. The original model also overpredicts the observed space density of CT AGN with 
logLx > 44 at 1.5 < z < 2, while underpredicting the z > 2.5 data points. The decadal 
CT AGN space density is shown in figure [2] where data and predictions for CT AGN with 
logLx = 44-45 are colored blue, and data and predictions for CT AGN with logLx = 45-46 
are colored green. The dashed lines again show the predictions of the original model. This 
model agrees fairly well with the logLx = 44-45 and logLx = 45-46 data at low redshift, 
but underpredicts the density of high redshift sources. As seen in figure [31 the original 
model agrees fairly well with the Swift/BAT log A^-logS" relation, where red lines denote the 
transmission dominated CT AGN number counts and dashed lines denote the original model. 
This model overpredicts the BAT number counts at the lowest flux levels but is in decent 
agreement with observations at the higher flux levels. A x 2 test was performed using the 
cumulative CT space density data points, the decadal CT space density data points which 
are not lower limits, and the BAT CT AGN number counts data points, giving 32 degrees 
of freedom. The original model was found to have a reduced x 2 °f X 2 re d = 15. 

We then attempted to model the CT population as only including sources with A > 
—0.05 independent of the fraction of type 2 AGN at any Lx and z. In order to match the 
peak of the CXB, ~91% of these AGN at any time must be CT. As seen in figures [TJ 121 and 
[3j where the predictions from this model are shown as the dot-dashed lines, this scenario 
underpredicted the observed local CT AGN space densities, overpredicted the Swift/BAT 
logAMogS' relation by a factor of 2.0 to 3.5, and resulted in x 2 e <i = 25. Thus, the CT 
population cannot be modeled as only including sources with very high Eddington ratios. 

The scenario where all CT AGN have A > —0.05 failed t o agree w ith local observations, 



and as local AGN are known to have low Eddington ratios (jHol 120081 ). low Eddington ratio 
sources were then added to the CT AGN population model. In order to fit the peak of 
the CXB, while not overestimating the CXB in the 100-200 keV range, this composite 
model predicts ~60% of AGN with A < -2.0 and ~86% of AGN with A > -0.05, at any 
time and independent of the fraction of type 2 AGN, are CT. The CXB fit found using 
this composite model is shown in figure H] with the contribution made by AGN accreting at 
different Eddington ratios plotted as separate lines. In this model the CXB is still dominated 
by obscured sources accreting at moderate Eddington ratios, i.e. type 2 Seyferts. The 
composite model predictions are shown as the colored solid lines in figures [U [21 and |3j 
This model is in fairly good agreement with the observed space densities for CT AGN with 
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\ogLx > 43 and logLx > 44, except for overpredicting the observed space density of CT 
AGN with logLx > 44 at 1.5 < z < 2. The composite model predictions match the 
observed cumulative space densities better than the original model predictions. At z < 1, 
the composite model slightly overpredicts the space density of CT AGN with log Lx = 44-45 
and underpredicts the observed densities at z > 2.5. For CT AGN with logLx = 45-46, 
this scenario agrees fairly well with observations. In general, the composite model is a better 
fit to the high luminosity, high redshift data than the original model. This scenario does 
overpredict the BAT number counts at the lowest flux levels, but is in decent agreement 
with observations at the higher flux levels. We find Xred = 9.8, indicating that the composite 
model is in better overall agreement with the observed CT AGN space densities and the 
Swift/BAT log iV-log S relation when compared to the original model, or one where all CT 
AGN have A > -0.05. 

Other models were investigated in which the CT population was assumed to be accreting 
at only low or moderate Eddington ratios. When the CT population was assumed to all 
be accreting at A < —2.0 the model underpredicted the observed high luminosity CT space 
densities and slightly underpredicted the Swift/BAT log N-\ogS relation. When —2.0 < A < 
—0.5 accretion was assumed for the CT population the model overpredicted the observed 
local space density by a factor of 3.2 and significantly overpredicted the Swift/BAT \ogN- 
logS relation. In conclusion, it is necessary to include both sources with A < —2.0 and 
A > -0.05 to model the CT AGN population. 



4. Discussion and Summary 



We find that current data suggests that the majority of AGN accreting at close to 
the Eddington limit are CT, in agreement with galaxy merger simulations (e.g., Hopkins 
et al. 2006) that show that mergers funnel gas into the central regions of galaxies igniting 
both starbursts and high Eddington ratio accretion onto the central black hole. These AGN 
are able to accrete rapidly because of a large abundance of gas and dust in the central 
region of the host bulge. This abundan ce of gas and dust naturally leads to a very high 
obscuring column density ( iFabianl Il999l ). Furthermore, radiat ion pressure from a rapidly 
accreting AGN will blow out any dusty gas with logA^H < 24 (IFabian et al. 1 120091 ) . Thus, 
in accordance with our findings, the vast majority of rapidly accreting, obscured AGN must 
be CT. 



A new result is that the CT AGN population, in addi t ion to including sources with 
A > —0.05, also includes sources with A < —2.0. While iHol (|2008[) explains that mos t ob- 
served low Eddington ratio sources tend to be unobscured. iTerashima fc Wilson ( 2003 ) find 
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evidence for a population of highly obscured low Eddington ratio AGN. Indeed iMalizia et al 
( 120 101 ) find an obscuration-luminosity relation in their sample of INTEGRAL sources which 
claims that lower luminosity sources tend to be more obscured. As these observed lower 
luminosity sources are local, these AGN are likely to have low Eddington ratios. Since low 
Eddington ratio AGN are only weakly accreting, they will have very little effect on their envi- 
ronment. Therefore molecular cl ouds would be able to come deep into the core of the bulge 
without being affected. Indeed iTerashima fc Wilson! (120031 ) find that the highly obscured 
low Eddington ratio sources in their sample might exhibit variable obscuration, as would 
be expected if the obscuration were due to molecular clouds near the black hole. Several 
studies have been conducted on so called 'changing-look' AGN (see Bianchi et al. 2005 and 
references therein), which are AGN whose spectrum changes from Compton thin to reflec- 
tion dominated on the time scale of years. It is believed that the changes in the spectrum 
of these local, low luminosity AGN is due to either the nucleus being obstructed by a cloud 
with a CT column density or the nucleus switching off, so that the re fl ection dominated 
state is actually the echo of a previous accretion episode. iBianchi et al. I (120051 ) find that it 
is most likely that the four AG N in their sample are in CT states. In a study of 82 LINERs, 
Gonzalez-Martin et al. I ( 120091 ) find about half of the LINERs in their sample appear to be 
CT and that there is a higher percentage of CT LINER s than CT Seyferts. The avera ge 
Eddington ratio of the CT LINERs in the sample used by iGonzalez-Martm et al. I ( 120091 ) is 
A = -4.7. This finding is consistent with the composite model for CT AGN discussed here 
which finds that ~60% of AGN with A < -2.0 are CT. 

It is clear that previous attempts to model the elusive CT AGN population as an 
extension of the Compton thin type 2 population fai l to exp lain the observed CT AGN 
space density and log iV-log S relation. iTreister et al. I ( 2009bJ ) suggest that CT AGN only 
contribute ~10% of the CXB and that the normalization of the CXB is overestimated by 
population synthesis models which predict a large population of CT AGN. However, changing 
the normalizatio n of the CXB does not account for the strong evolution of CT AGN at z = 
1.5-2.5 found by ITreister et al. I (l2009af). Here we assume that CT AGN contribute ~20% 
of the CXB based on the findings of iDraper fe Ballantynd ( 120091 ). but the composite model 
presented here does account for the strong evolution of CT AGN at high redshift while not 
overestimating the local CT AGN space density. Additionally, models which include CT 
sources with —2.0 < A < —0.5 overpredict locally observed CT AGN space density. Models 
which only include CT sources with A < —2.0 underpredict the density of high luminosity 
CT sources. If CT sources are assumed to only have A > —0.5, the model overpredicts the 
observed log iV-log S relation. A composite model, where the CT AGN population includes 
both sources accreting at close to Eddington and sources accreting at <1% of Eddington, 
best explains observations. 
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As CT AGN are rare and difficult to observe due to high obscuration, it is difficult to use 
observations to directly decipher the physical phenomenon which gives rise to the extreme 
levels of obscuration found in CT AGN (e.g., Triester et al. 2009b; Rigby et al. 2009). 
Thus physical models of CT AGN must use indirect observational constraints, like observed 
space densities and number counts, to uncover the nature of the extreme obscuration of 
these sources. This letter presents a model which constrains the physical parameter of the 
Eddington ratios of CT AGN to >90% and <1%. In the future the understanding of CT 
AGN will be expanded through IR and sub-mm studies of the reprocessed radiation from 
CT AGN and very hard X-ray imaging with missions like NuSTAR. 

The authors thank A. Merloni for helpful conversations. 
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Fig. 1. — Cumulative CT AGN space density. Red points and lines denote CT AGN with 
logLx > 43 and blue points and lines denote CT AGN with logLx > 44. Solid lines show 
the composite model where CT AGN have either A < —2.0 or A > —0.05. The dashed lines 
show the original model where CT AGN are considered a simple extension of the Compton 
thin type 2 population. The dot-dashed lines show the scenario where all CT AGN have 
A > —0.05. Data points are shown from several studies: triangles: Treister et al. (2009a); 
squares: Tozzi et al. (2006); stars: Alexander et al. (2008); pentag ons: Fiore et al . (2009) ; 



circle: local Swift/BAT and INTEGRAL data point reported by Treister et al. ( 2009bl) 



adjusted to reflect the flux-luminosity relation for CT AGN described by lRigby et al. I ( 120091 ). 
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Fig. 2. — Decadal CT AGN space density. Blue points and lines denote CT AGN with 
log Lx = 44-45 and green points and lines denote CT AGN with log Lx = 45-46. Lines are 
the same style as in figure [TJ Data points are shown from several studies: filled squares: Yan 
et al. (2007); circles: Martinez-Sansigre et al. (2006); stars: Polletta et al. (2006); triangles: 
Fiore et al. (2009); pentagons: Fiore et al. (2008); open squares: Alexander et al. (2008). 
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Fig. 3. — BAT number counts. Blue lines signify total AGN number counts and red lines 
signify transmission dominated CT AGN number counts. Colored lines are the same style 
as in figure HJ The black line shows the total AGN log iV-log S relation from Tueller et al. 
(2008) and the black data points show the CT log A-log S data points from Treister et al. 
(2009b). 
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Fig. 4. — CXB fit using the composite model showing contributions of AGN accreting at 
different Eddington ratios. The solid black line shows the combined contribution of non- 
blazar AGN and blazars to the CXB. Colored lines show the contribution of AGN accreting 
at the color designated Eddington ratio, where dashed lines refer only to the Compton thin 
contribution and solid lines show the combined contribution of Compton thin AGN and CT 
AGN. Measurements from various instruments are shown as colored areas and data points; 
blue: ASCA GIS (Kushino et al. 2002); magenta: Rossi X-ray Timing Explorer (RXTE; 
Revnivtsev et al. 2003); green: XMM- Newton (Lumb et al. 2002); red: BeppoSAX (Vecchi 
et al. 1999); yellow: ASCA SIS (Gendreau et al. 1995); cyan: XMM- Newton (De Luca & 
Molendi 2004); grey data: HEAO-l (Gruber et al. 1999); blue data: INTEGRAL (Churazov 
et al. 2007); red data: Swift/B AT (Ajello et al. 2008); black data: Swift/XRT (Moretti et 
al. 2009); green data: INTEGRAL (Tiirler et al. 2010). 



